Is the shell-focusing singularity of Szekeres space-time visible? 
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The visibility of the shell-focusing singularity in Szekeres space-time - which represents quasi- 
spherical dust collapse - has been studied on numerous occasions in the context of the cosmic 
censorship conjecture. The various results derived have assumed that there exist radial null geodesies 
in the space-time. We show that such geodesies do not exist in general, and so previous results on 
the visibility of the singularity are not generally valid. More precisely, we show that the existence of 
a radial geodesic in Szekeres space-time implies that the space-time is axially symmetric, with the 
geodesic along the polar direction (i.e. along the axis of symmetry). If there is a second non-parallel 
radial geodesic, then the space-time is spherically symmetric, and so is a Lemaitre-Tolman-Bondi 
(LTB) space-time. For the case of the polar geodesic in an axially symmetric Szekeres space-time, 
we give conditions on the free functions (i.e. initial data) of the space-time which lead to visibility 
of the singularity along this direction. Likewise, we give a sufficient condition for censorship of 
the singularity. We point out the complications involved in addressing the question of visibility of 
the singularity both for non-radial null geodesies in the axially symmetric case and in the general 
(non-axially symmetric) case, and suggest a possible approach. 
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I. INTRODUCTION AND SUMMARY. 

The weak cosmic censorship hypothesis (CCH) main- 
tains that rcahstic gravitational collapse leads to the for- 
mation of a black hole rather than a naked singularity. 
Among the different studies of the hypothesis, we men- 
tion two categories of interest. First, there are demon- 
strations of the validity of the hypothesis in specific cir- 
cumstances (the prime example of this is Christodoulou's 
proof of the instability - and hence non-realistic nature - 
of naked singularities in the spherical collapse of a min- 
imally coupled scalar field The other category in- 
volves the construction of an example of a space-time 
which undergoes collapse from a regular configuration to 
a naked singularity. Many such examples have been con- 
structed, but to date, none has been shown to involve 
both (i) a physically realistic matter model and (ii) sta- 
bility in the initial data space of those space-times which 
give rise to naked singularities. It is probably fair to say 
that the main utility of the latter class of studies has 
been to refine and better understand the content of the 
CCH. 

These examples have mainly involved spherically sym- 
metric space-times, for example the shell-crossing [2, Q 
and shell-focusing [3| singularities in Lemaitre-Tolman- 
Bondi (LTB) spherical dust collapse, and the naked sin- 
gularity solutions that arise at the threshold of black hole 
formation in scalar field and perfect fluid collapse Q. A 
notable (although as we will see flawed) exception to this 



involves the various studies of the visibility of the shell- 
focusing singularity in Szekeres space-time [6j-[l^. First 
analyzed by Szekeres in [Tsj , this class of space-times cor- 
responds to solutions of the Einstein equations for dust, 
where the fluid flow vector is geodesic and non-rotating. 
The metric admits no Killing vector fields (in the general 
case) but for reasons described below is referred to as 
quasi-spherical. It can be understood as a non-spherical 
generalization of the LTB class of dust-filled space-times, 
and its evolutionary aspects are very closely related to 
those of the corresponding LTB models, and so are rel- 
atively straightforward. Hence analyzing the visibility 
or otherwise of singularities that arise in this model af- 
fords an opportunity to study cosmic censorship in non- 
spherical collapse. 

Szekeres initiated the study of the singularities in this 
model, and noted the possibility that the shell-crossing 
singularity (see below) may be visible 0|. He also 
noted the occurrence of an apparent horizon (or in cur- 
rent terminology, a marginally trapped tube) that forms 
at least as early as the shell-focusing singularity (again, 
see below). As in the spherical case, the shell-crossing 
singularity is interpreted as being fundamentally non- 
gravitational in origin and not of particular significance 
for cosmic censorship. Thus attention turned to the shell- 
focusing singularity, where in the spherical case, a more- 
or-less complete understanding of the visibility or other- 
wise of the singularity has been developed. In these stud- 
ies @-|l2|, the question of visibility of the shell- focusing 
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singularity is considered from the point of view of the 
existence or otherwise of future-pointing null geodesies 
with an additional simplifying property that allows one 
to refer to these geodesies as 'radial' (see below). We 
show below that radial geodesies do not exist in a gen- 
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eral Szekeres space-time and consequently, the anal yses 
of the question of visibility of the singularity of fdil-fl^l 
are not generally valid. Motivated by this observation, we 
revisit the question of the visibihty of the shell-focusing 
singularity. We provide some preliminary results on this 
question, and in particular, consider it in the case when 
the space-time is axially symmetric. In this case, a single 
radial geodesic direction exists, and the analysis of the 
visibility of the singularity along this direction is essen- 
tially the same as the spherically symmetric case. We 
emphasize that the question becomes considerably more 
difficult in the general (non-axially symmetric) case. 

In the next section, we review the basic properties of 
Szekeres space-time, and discuss what is meant by refer- 
ring to this space-time as 'quasi-spherical'. We then dis- 
cuss the formation of singularities in a collapsing Szekeres 
space-time, and discuss the conditions on the free met- 
ric functions that arise from the imposition of regularity 
conditions on the initial data. Here, we specialize to the 
so-called marginally bound case. We will indicate clearly 
when this restriction is in place, and when results apply 
generally. 

In section 3, we analyze the geodesic equations of Szek- 
eres space-time and show that the existence of a radial 
geodesic implies that the space-time is axially symmet- 
ric. Furthermore, we show that the existence of a second 
non-coUinear radial geodesic implies that the space-time 
is spherically symmetric. The marginally bound assump- 
tion is not required in this section. 

In section 4, we derive some elementary results in the 
marginally bound case relating to the visibility of the 
shell-focusing singularity, which is spherically symmetric 
in the sense that it corresponds to a surface t = tc{r) , r > 
0. This singularity is always preceded by an apparent 
horizon given hy t ~ tah{i^)'- we have tahij") < ic(?') for 
all r with equality if and only if r = 0. The region of 
space-time with tahi^) < t < tc{r) is trapped. Then 
intuitively, one expects that only the central singular- 
ity {t,r) = (<c(0),0) can be visible. This is immediate 
in spherical symmetry, but is slightly non-trivial in the 
quasi-spherical case - the result is confirmed nonetheless. 
We derive the related result that a geodesic that emerges 
into the future from the central singularity must emerge 
into the untrapped region t < tah, and hence show that 
the singularity is censored if t'^f^{r) < (this condition 
can be described in terms of the initial data of the col- 
lapse). We also show that for sufficiently small values 
of r, the apparent horizon is a one-way membrane for 
all future-pointing causal geodesies: such geodesies can- 
not leave the trapped region. Again, this is immediate 
in the spherical case, but requires checking in the quasi- 
spherical case. 

In section 5, again working in the marginally bound 
case, we consider the question of the visibility of the sin- 
gularity for the polar null geodesic of the axially sym- 
metric models. This problem is essentially the same as 
the corresponding spherically symmetric problem, and 
we can give sufficient conditions (in terms of the initial 



data) for the formation of a naked singularity. 

We conclude by briefly considering the substantive, 
and crucially, open, question of the visibility of the shell- 
focusing singularity in quasi-spherical collapse. We point 
out how this question is considerably more difficult than 
in the spherical case and suggest an approach to its con- 
sideration. We set 8ttG = c = 1. 



II. SZEKERES SPACE-TIME AND ITS 
SINGULARITIES. 

A comprehensive review of the properties of the Szek- 
eres space-times representing non-accelerating, irrota- 
tional dust is given in [Tsj : this review includes equiva- 
lent invariant characterizations of the class which involve 
some technicalities that will not play any role here. Suf- 
fice to say that these lead uniquely to the line element 



dt^ + e^"dr^ + e^f^idx^ + dy^) 



(1) 



where 



= R{t,r)e'', (2) 
e-'' = A{r){x^ +y^) + Bi{r)x + B2{r)y + C{r), (3) 
R' + Rv' 



(4) 



where the prime denotes diff'erentiation with respect to r. 

d 

The geodesic fluid flow vector is — and the coordinates 

at 

(r, x, y) are co-moving. The ranges of these coordinates 
are (a;,y) € K^, r > and that of t will be discussed 
below. The free fmictions A, Bi, B2 and C are related by 



Bf + B^- AAC = -1. 



(5) 



(There are also solutions with and +1 on the right hand 
side here: the choice — 1 forms part of the input essential 
to the interpretation of these space-times as being quasi- 
spherical.) The remaining Einstein equations determine 
the evolution of R and define the energy density of the 
space-time: 



F{r) 
R 



p(t, r, X, y) = 



F' + iFv' 
R^{R' + Ru' 



(6) 



(7) 



F > is a function of integration. We note that 
the coordinate freedom corresponds to rcscalings of the 
co-moving radial coordinate r — > f{r) (which must 
be a monotone mapping), the transformations in the 
(x,?/)— plane discussed in ([8]) below and trivial shifts of 
the origin of t. 

The quasi-spherical interpretation arises as follows 
II4I : each 2-surface St r of constant t and r is a round 
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2-sphere with proper radius R(t,r). To see this, we write 
the hnc element dl'^^ ^.j of St^r in the stereographic coor- 
dinates C, — X + iy: 



dl 



R\t,r) 



(Air)a + B{r)C + B{r)C + Cir)y 



rdCdC 



where B ~ {Bi — iB2)/1. The form of this hne element 
is invariant under the fractional linear transformation 



c->e = 



mC, + n ' 



kn — Im = 1 



(8) 



and for each fixed value of t and r, such a transformation 
can be found so that the line element has the form 



which is the line element of the round 2-sphere with ra- 
dius R. It should be noted that the condition ^ plays 
a crucial role in deriving the explicitly spherical form of 
d^ft r)- transformation (x, y) (6, (jj) given by 

C = e^^cot^ 

yields the more familiar spherical form 

= R^{t,r){dd^ + sin^ 0d(j)^). 

The fact that a different transformation ([8|) is required 
for each different St^r indicates that in each 3-space of 
constant t, these 2-spheres are not concentric. 

It will be useful to consider the form of the 4- 
dimensional line element using the spherical coordinates 
(9,(1)). This yields 

ds^ ^ -df + e^"'dr'^ + e^^* {dd^ + sin^ edct)"^ ), (9) 

where now 



e'^ = R{t,r)e^ 



(10) 



a cos 6* + 6i sin6'cos(/) + &2 sin 6' sin + c, (11) 



R' + Rfi' 



(12) 



VT+7 

The functions a, hi, c are related to A, Bi,C by 

a^A-C, h^B,, c = A + C. (13) 

The condition ([5]) reads 

-a"^ -bl-bl = 1. (14) 

Before proceeding to discuss gravitational collapse and 
the formation and nature of singularities in these space- 
times, we note that the spherical limit arises when and 
only when Bi = B2 = and A and C are constant and 
equal (and so by ([5]) both equal to 1/2). Furthermore, it 



is clear from Q and (llip that the space-time is axially 

d 

symmetric with Killing vector — when Bi = B2 = 0. 

In order to model gravitational collapse using a Szek- 
eres space-time, we choose the negative root of ©. The 
resulting equation, and the analysis of its consequences, 
are greatly simplified by taking /(r) = 0. By anal- 
ogy with the spherical case, this is referred to as the 
marginally bound case. The equation is then easily inte- 
grated and the solution can be written in the form 



(15) 



where tc is a function of integration that describes the 
time at which the 2-sphere St^r collapses to zero radius. 
This is called the shell-focusing singularity: all the 'shells' 
St^r collapse to zero radius at this surface. 

It is convenient to exploit the freedom in the co-moving 
radial coordinate r to set i? = r on an initial surface: by 
an allowed shift of the origin of t, we can take this to be 
t = without loss of generality. Thus R{0, r) ~ r and 
hence i?'(0,r) — 1. With this choice, we have 
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ic(r) = -^-, 
which leads to the following convenient form of ([15 



(16) 



(17) 



It can be shown that the Kretschmann scalar of ^ 
diverges if and only if the density ([7]) does so. Thus scalar 
curvature singularities can be discussed in terms of this 
quantity alone. As well as the shell-focusing singularity 
at i = ^c(^) for which i? = 0, there may be singularities 
when R' + Rv' vanishes. By analogy with the spherical 
case, these are referred to as shell-crossing singularities. 
(The analogy is perhaps not quite appropriate: shells of 
Szekeres space-time, i.e. the 2-spheres St.r will cross if we 
encounter i?(t, ri) = i?(t, 72) for some ri 7^ r2 and some 
t > 0. A necessary and sufficient condition for this to 
occur is R'{t, r) = at some t > 0, r > 0.) A crucial task 
is to rule out the occurrence of both types of singularity 
on the initial slice. To this end, we note that the initial 
density is given by 

po r := P(0> n = -— — -. 

r^(l + rv') 

We will require this term to be non-negative and finite 
for all r > 0. Thus we impose 



F' + 2,Fu' > 0, 



1 + riy' >0 



(18) 



(19) 



for all r > and all {x,y) e R^. We note that we also 
impose e"" > 0. Identical conditions hold with 1/ re- 
placed by fi. These conditions - i.e. ([TO)) and > 
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- have been considered by Szekeres [ij], and it is worth 
repeating the result here as we will use the corresponding 
conditions on a,bi and c below: 

Lemma 1 (i) Let the condition (0j hold. Then e^^ = 
Aix^ + j/^) + Bia; + i?2y + C is positive for all {x, y) <E M-^ 
and all r > if and only if 

A{r) > for all r > 0. (20) 

(a) Assuming the conditions 0) and e^"^ > 0, 

1 + ru' > 

for all {x,y) E M.^ and all r > if and only if 

A-rA' >Q for all r > (21) 

and 

A'C - B'B' >-\ for all r > (22) 



Requiring that pq be finite in the limit as r — s- 0, Szek- 
eres also derives the condition 



F{r) = 0(r 







(23) 



space-time will contain shell-crossing singularities in its 
evolution. 

To summarize, we assume the conditions on A, Bi and 
C of Lemma 1 and the conditions and ([M]) on F. 
These guarantee that the collapse proceeds from a regular 
state to the shell-focusing singularity, and that no shell- 
crossing singularities occur prior to the shell-focusing sin- 
gularity. 



III. GEODESICS, RADIAL GEODESICS AND 
SYMMETRY 



With a clear description of the shcU-focusing singular- 
ity in place, we can now consider its causal nature or 
more accurately, the question of its visibility. Of course 
this requires the analysis of geodesies, and the determi- 
nation of whether or not there are future pointing causal 
geodesies that emerge from the singularity. We will use 
a subscript {at etc.) to denote partial derivatives with 
respect to t, x and y. We retain the prime for derivatives 
with respect to r, and an overdot will represent differen- 
tiation along the geodesic (e.g. with respect to an affine 
parameter for null geodesies). We also note that we can 
drop the assumption that the space-time is marginally 
bound. The geodesic equations for the line element ([T]) 



which we shall assume henceforth. 

Next, we will consider conditions that rule out the oc- 
currence of a shell-crossing singularity [R' -I- Rv' ~ 0) 
prior to the occurrence of the shell-focusing singularity. 
That is, we seek conditions on the initial data functions 
so that 

R' + Rv' > 0, for all 0<t< t^ir), r > 0. 



We note that this condition follows by (|T9l) if the in- 
equality R' — R/r > holds. But using (fTT]) . this latt er 
condition is equivalent to > 0. Hence the condition 



F 

F'-3- <0, 
r 



r > 0, 



(24) 



which is equivalent to > 0, is a condition that can 
be imposed on the initial data and that guarantees the 
absence of shell-crossing singularities. It is worth noting 
that for i > 0, p4|l is equivalent to the condition 



R'>^. 

r 



(25) 



In a sense, this is the best bound that can be imposed. 
First, it includes all cases of interest: tj, > is a necessary 
condition for the visibility of the shell-focusing singular- 
ity. (To see this, we simply note that t must increase 
along a future pointing causal geodesic emerging from the 
singularity.) Second, if the bound ([M|) is violated, then 
we can find examples of v (which is constructed entirely 
from initial data functions) for which the corresponding 



i+ate^'^r^ + Pte^'^ix^+f) = 0,(26) 
f -|- a'r^ + 2atir + 2axXr + 2ayyr 

-/?'e2/3-2"(i2^y2) ^ 0,(27) 



X + Po^ii'^ - f ) + '^Pyiy + '2l3tix + 2(3' rx 
y + Pyif - i') + 2/3.iy + 2Ptiy + 2(5' fij 



0, (28) 



„2q-2/3-2 _ 



= 0, (29) 



and we have the first integral 



(30) 



where e = 0, +1,— 1 for null, space-like and time-like 
geodesies respectively. 

In radial geodesies are defined to be those 

along which x and y have constant values. From (pS]) 
and ([29]) . we see that along such a geodesic, we must 
have 

^g2a-2;3 -2 ^ ^g2o-20 .2 ^ q (3^) 

OX ay 

We rule out e°~'' = 0, as this corresponds to a singularity 
(and the geodesies must reside in the space-time rather 
than on its singular boundary). Both equations of (|3ip 
are satisfied if we take f = 0. The only possible solutions 
of the geodesic equations then have P = 1 = — e: these 
are the fiuid fiow lines and we will refer to these as trivial 
radial geodesies. In particular, there are no ttoH geodesies 
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that satisfy the condition f = 0. The only possibihty that 
remains in (j3ip is that 

da da ^ 
dx dy 

along the geodesic. It is clear that this is a restriction 
on the metric functions. We can determine the exact 
geometric nature of this restriction. 
In order to do so, we note first that 

da _ Re-"' d'^v 
dx \/l~+T drdx' 

with a similar result holding for ay. Then a straightfor- 
ward calculation shows that the vanishing of is equiv- 
alent to the vanishing of 

Qi := {AB[-A'Bi){x^-y^) + 2{AB!,-A'B2)xy 
+2{AC' - A'C)x + {BiB'2 - B[B2)y 
+ {BiC'-B[C), (32) 

while vanishing of ay is equivalent to vanishing of 

Q2 := [AB'^- A'B2){y'' -x^) + 2{AB[- A'Bi)xy 
+2{AC' - A'C)y + {B2B[ - B'^Bi)x 
+ {B2C' - B'^C). (33) 

We can now prove the following result. 

Proposition 1 If a Szekeres space-time admits a non- 
trivial radial geodesic, then it also admits a Killing vector 
field generating an axial isometry. 

Proof: Suppose that there is a non-trivial radial geodesic 
along which {x,y) = (xo^yo) is constant. Using a trans- 
formation of the form ([5]), 

(, = X ty ^ t, = u + tv = — 

mC + ri 

we can assume without loss of generality that (xo,yo) ~ 
(0,0). The necessary conditions Qi = Q2 = for the 
existence of a non-trivial radial geodesic then yield 

Bi = AiC, B2 = A2C 

for some constants Ai, A2. (We note that ^ implies that 
C ^0.) If we consider a further coordinate transforma- 
tion of the form ([8]) - but with / = to preserve the origin 
- we find that (with B = [Bi - iB2)/2)) 

B ^ B^ = kfiB + mnC = n{kX -t- m)C 

where A ~ (Ai — iX2)/2. Since kn — Im ~ kn ~ I, we 
have k ^ ^ n, and so we can choose m = ~kX to get 
= 0. Thus by using the coordinate freedom in the 
stereographic coordinates {x,y), we can write the line 
element ([T]) in a form in which Bi = B2 — 0- As seen 
in Section 2 above, this is a sufficient condition for the 
space-time to be axially symmetric, with axial Killing 
field given by ^ = -y^ +x-§^. ■ 



Remark 1 In this Proposition, the radial geodesic is 
normal to each of the 2-spheres and emerges from the 
point with stereographic coordinates x ^ y = 0. In 
spherical coordinates, this corresponds to 9 = n: the 
south pole in the standard configuration. The north pole 
(6 ~ 0) corresponds to the point at infinity in stereo- 
graphic coordinates, and so is not covered by the coor- 
dinate patch {x,y) S M^. However this point can be 
included by using an additional coordinate patch and it 
is then clear that there is also a radial geodesic emerg- 
ing from the north pole. We will refer the these coUinear 
radial geodesies as the polar geodesies. 

Corollary 1 If a Szekeres space-time admits two non- 
collinear non-trivial radial geodesies, then the space-time 
is spherically symmetric. 

Proof: From the previous proposition, we may situate 
the first non-trivial geodesic in the direction (x, y) = 
(0,0). Then as we have seen Bi = B2 = and the 
space-time is axially symmetric. A second non-collinear 
non-trivial radial geodesic has {x,y) = {xo,yo) 7^ (0,0) 
constant along the geodesic, and as we have seen, a neces- 
sary condition for the existence of such a geodesic is that 
Qi and Q2 vanish along the geodesic. As i?i = i?2 = 
and xq and j/q are not both zero, (|32p and ([33]) yield 

AC - A'C = 0. 

The condition ([5]) in the present case gives AC = 1/4. 
Combining this with the previous relation shows that A 
and C are both constant. A transformation of the form 
([5]) (with ^ = to preserve the origin) can then be used 
to set 2 A ~ 2C = 1, and so the line element is spherically 
symmetric. ■ 

Remark 2 We can give a geometric interpretation of 
this corollary. Proposition 1 shows that to each non- 
trivial radial geodesic of Szekeres there corresponds an 
axis of symmetry. Then the existence of a second 
non-collinear non-trivial radial geodesic implies that the 
space-time admits two non-parallel axes of symmetry and 
so must be spherically symmetric. 

Remark 3 The results of this section imply that the 
analysis of which were carried out for radial null 

geodesies, can only be valid for the polar geodesic of an 
axially symmetric Szekeres space-time, or for a spher- 
ically symmetric Szekeres space-time - i.e. LTB space- 
time. Thus the question of the visibility of the shell- 
focusing singularity in a general Szekeres space-time re- 
mains open. It is worth noting that the imposition of 
the assumption of axial symmetry, equivalent to setting 
Bi ~ B2, implies that the number of free functions 
{A, Bi,B2,f, F) has been reduced from 5 to 3. Therefore 
this sector of the Szekeres class is highly specialized - one 
could impose a topology on the space of free functions in 
which the non-axially symmetric solutions comprise an 
open dense subset of the whole space - and so the results 



6 



pertaining to the polar geodesies cannot be assumed to 
reflect the general behaviour. It is possible that the naked 
singularities found in the axially symmetric case along 
the polar direction are (i) not visible from any other di- 
rection and/or (ii) are not present in the non-symmetric 
case. 



which is positive by the no-shell crossing singularity con- 
dition.) As all stationary points must be local minima, 
there can in fact be only one local minimum. Thus 
/■(t) < for all T < tq. Suppose that this geodesic were 
to meet the singularity t = tc{r). Now along the geodesic 
we have 



IV. SOME BASIC RESULTS IN THE 
MARGINALLY BOUND CASE. 

We have argued above that the question of cosmic cen- 
sorship in Szekeres is open. In the remainder of this pa- 
per, we seek to address this question. In this section, we 
will derive some results, valid in the marginally bound 
case only, which provide some useful preliminary results 
for the study of the visibility of the shell-focusing singu- 
larity. These results share the feature that they are trivial 
in the spherically symmetric case: in the quasi-spherical 
case, some checking is required. 

To begin, we recall from the work of Szekeres 
that an apparent horizon forms at the hypersurface 
R{t,r) = F{r). That is, the outgoing future-pointing 
null geodesic normals to each 2-sphere St,r have zero ex- 
pansion on this hypersurface. By outgoing, we mean that 
r increases with the afRne parameter along the geodesic. 
Note also that while these geodesies are initially and in- 
stantaneously radial, this condition is immediately vio- 
lated as the geodesic moves away from the 2-sphere to 
which is was normal - see ([28]) and ([29|) . 

From (fT7|) . we can show that the apparent horizon is 
given by 



Rti + R'f 



t = tah{r) = tc{r) 



iFir). 



(34) 



Thus the apparent horizon precedes the shell-focusing 
singularity (F > for ?- > 0), and the condition (IM)) im- 
plies that the apparent horizon and the shell-focusing sin- 
gularity meet at the central singularity {t, r) ~ (ic(0): 0). 
As shown in the region tahir) < t < tc{r) , r > is 
trapped. That is, the 2-spheres St^r are closed trapped 
surfaces in this region. The region t < tah{r),r > is 
untrapped. 

First, we point out that the portion r > of the shell- 
focusing singularity is censored. 

Proposition 2 There are no future-pointing causal 
geodesies of a marginally bound Szekeres space-time with 
past endpoint on the surface t = tc{r) for any r > 0. 

Proof: We consider a future-pointing {t > 0) outgoing 
(r > 0) causal geodesic of Either the geodesic re- 
mains outgoing, or we encounter a value Tq of the geodesic 
parameter r (affine parameter or proper time) for which 
r(To) = 0. But then (|27p shows that this is a local min- 
imum (f > 0) of r along the geodesic. (To see this, we 
note that 



dr 




(35) 



Since t > 0, R' > and r < 0, this must by negative 
in the approach to the singularity. But the singularity 
occurs at i? = 0, and so R cannot increase into the past 
{R < 0) to reach the singularity, and so we get a contra- 
diction. 

We can now assume that r > for all t < tq where tq 
is some arbitrary initial value for r. This being the case, 
we can use r as a parameter along the geodesic. Then 
along the geodesic we may use ([30| to write 



dt 
dr 




(36) 

The fact that the geodesic is both future pointing and 
outgoing indicates that the correct (positive) root has 
been taken here. The change of R along a future pointing 
geodesic that emerges from the shell-focusing singularity 
at some r > satisfies the following: 



F dt , 
RTr^^ 




R' = -Rv' 



(37) 



(i?' + Rv')e'' 



The second line comes from the field equation Q in the 
marginally bound case. The third line arises due to the 
fact that for sufficiently small R. the geodesic must be 
in the trapped region for which F > R. The fourth line 
follows from p6|) above and the definition ([4]). The last 
line comes from the initial regularity condition (jl9p . In- 
tegrating the overall inequality proves the stated result: 
R cannot reach zero (its value on the shell-focusing sin- 
gularity) unless r also drops to zero. ■ 

Corollary 2 A future pointing causal geodesic with past 
endpoint on the shell-focusing singularity must have its 
past endpoint on the central singularity (t, r) = (tc(0)) 0)- 



Corollary 3 A future pointing null geodesic with past 
endpoint on the central singularity must emerge into the 
untrapped region of space-time. 
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Proof: Suppose on the contrary that the geodesic 
emerges into the trapped region, i.e. there exists 6 > 
such that R{t{r),r) < F{r) for values of R along the 
geodesic and for all Q < r < S. Letting tq € (0, S) and in- 
tegrating (|37p from r to vq shows that in the R-r plane, 
R stays above the line R = for all r < rg. Using 
this implies that i? > F for sufficiently small values 
of r, yielding a contradiction. ■ 

Proposition 3 If t'^fi{r) < on [0,6) for some 6 > 0, 
then the central singularity is censored. 

Proof: The proof follows immediately from Corollary 3: 
if a geodesic were to emerge from the central singularity, 
then it must emerge into the untrapped region and must 
have i'(r) non-negative for all sufficiently small values of 
r. This cannot happen if i',^ is negative in a neighbour- 
hood of the singularity. I 

Remark 4 We note that this repeats in the Szekeres 
case a result that holds in some generality in spherical 
symmetry [l6j . and provides a sufficient condition, in 
terms of initial data, for the singularity to be censored. 

Finally in this section, we prove another result that 
mirrors precisely the situation in the spherical case. As 
in that case, the proof relies crucially on some of our 
assumptions about the regularity of the initial data. 

Proposition 4 For sufficiently small values ofr, the ap- 
parent horizon t = tahif) o,cts as a one way membrane: 
a future pointing null geodesic cannot cross the horizon 
from the trapped to the untrapped region. 

Proof: Let p be a point of space-time on the apparent 
horizon with r > and consider a future-pointing null 
geodesic at p. If r|p = 0, then the fact that 

is finite for all r > and that t\p > proves the stated 
result. Suppose then that r\p ^ 0. Then there is a neigh- 
bourhood / 9 So such that f(s) ^ for all s G / with 
s\p — So where s is the parameter along the geodesic in 
question. Then for points on the geodesic correspond- 
ing to /, we can take r to be the parameter along the 
geodesic. Along a future pointing outgoing null geodesic, 
we then have, using 



Then using = ~ f-^'i show that 



# > 
dr 



> R'- 



R 



(38) 



where we have used (|4]) and the no-shell crossing condi- 
tion ((25|) . Using (fT7|) . we can show that 



t=ta.h{r} 



(39) 



r 



> t', ^ F' - - > 



The initial regularity condition (|23p indicates that this 
last inequality holds for sufficiently small values of r. 
Thus when projected onto the r-t plane, for sufficiently 
small values of r, an outgoing null geodesic can only cross 
the apparent horizon from below. The same result is im- 
mediate for ingoing null geodesies as we can show that 
subject to ([23|) . the slope of the apparent horizon is 
positive for small values of r. M 

Remark 5 We note that the result above is equivalent to 
stating that the apparent horizon is space-like for small 
values of r. Globally, there is no restriction: the hori- 
zon may be null or time-like for larger values of r, and 
can change character. Thus the apparent horizon is not 
always a one-way membrane in Szekeres space-time. 



V. POLAR GEODESICS IN THE AXIALLY 
SYMMETRIC CASE 

As we have seen in Section 3, the only case in which 
radial geodesies exist in Szekeres space-time is when the 
space-time is axially symmetric and that furthermore the 
only radial geodesies that can emerge are in the polar di- 
rection. In this situation, the analysis of the visibility 
of the singularity is essentially the same as that for the 
spherically symmetric case. We show here that there 
are choices of the initial data for which the central sin- 
gularity is visible along the polar direction. We follow 
the treatment of the spherically symmetric (LTB) case 
given in . However as with the previous section, some 
care must be taken to account for the minor differences 
between the present case and the spherically symmetric 
case. 

The axially symmetric case is obtained by setting Bi = 
B2 = in Q, and the polar geodesic corresponds to 
X = y = 0. Then the null geodesic equations (j26p - ([30|) 
reduce to 



t + att = 
f+{a' + 2ate")r'^ = 



(40) 
(41) 
(42) 



Our first step is to show that we can replace the affine 
parameter s by the coordinate r, and consider the projec- 
tion of the geodesic into the r-t plane. To see this, we note 
that from ((42|) . if r(so) = for some sq, then f(s) = for 
all s, and the geodesic reduces to a single point. So we 
can assume that r ^ along the geodesic. Hence a polar 
null geodesic that is initially outgoing r(so) > remains 
outgoing for all s. Consequently, apart from the question 
of the maximal s— interval of existence of the geodesic. 
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all information regarding the outgoing polar geodesic is 
contained in the single equation 



dt 
dr 



R' + Rv'. 



(43) 



Along such a geodesic 7, we have 

C 



V U = - 



C 



In the axially symmetric case, the conditions ([S]) and 
(Iini)-(I211) reduce to 



A>0, 
A-rA' > 0, 
1 



A'C > 



4^2 ■ 



From these we obtain 



- < ly U < -. 

r r 



(44) 

(45) 
(46) 

(47) 



(48) 



Thus the additional symmetry in the problem yields a 
useful additional bound on i^' {cf. (fT9|) ). 

To proceed, we make an additional mild assumption on 
the structure of the function F. We define the number 
/o and the function Fi by 

^^ = r3(/o + Fi(r)), Fi(0) = 0. 

Our mild assumption is that /o > 0: this corresponds to 
the initial central density being strictly positive. The no- 
shell crossing condition corresponds to F[ < for r > 0, 
and so /i := Fi{0) < 0. 

Proposition 5 If fi < 0, then there is a future pointing 
outgoing polar geodesic with past endpoint on the central 
singularity. 

Proof: For /i < 0, we can use ([34|) to write 

tah^^fa - g/o fir + 0{r ), 

where this and all other asymptotic relations in the 
present proof refer to the limit r —^ 0. For constant k, 
with < K < —hfn^^^ fi =: A, define 



Then for sufficiently small 5i , there is a non-empty region 

n[6ui^] = {{t, r) : t,{r) < t < tah{r), 0<r< 61}. 

Note that t'^ = k > 0. Along a future pointing outgoing 
polar geodesic 7, we have 

^=R' + R,y' <R' + -, 
dr r 



and a straightforward calculation using p?|) then yields 



dt 
dr 



< T 



1/3 



n2/3 



2/1 



,,2/3 



+0(r5/2). 



Since /i < and < k < A, the leading coefficient here is 
positive, and so there is a (52 > such that for all r < S2, 
we have 



dt 
dr 



< t' 



Hence by taking S to be sufficiently small to allow use 
of Proposition 4 and to minimize Si and 62^ we see that 
a future-pointing outgoing polar null geodesic in the re- 
gion ri[(5, k] as defined above cannot leave this region as 
we extend back into the past. Hence the geodesic must 
extend back to the central singularity r = 0,t = tc{0) — 

tahio) = t,io). m 

Remark 6 We note that as in the spherically symmet- 
ric case, it is possible to consider the case where /i = 0. 
This requires an additional assumption on the differen- 
tiability of F and on the value of the coefficients of a 
Taylor expansion of the function around r = 0. There 
is nothing to indicate that the results obtained in this 
way would differ from those obtained in the spherically 
symmetric case. 

We consider next the important question of whether 
or not these geodesies meet the singularity at some finite 
affine parameter value in the past, or if s — > — 00 as 
r — > 0,t — > tc{0). In order to do this, we study the 
sign of at in the region Cl[S, k] introduced in the proof of 
Proposition 5 above. 

Lemma 2 There are values of k G (0,— |/q ^^'^ fi) and 
(5 > such that at > Q for all (t, r) G fl[d, k\. 

Proof: We have 

R't + Rti^' 



at 



R' + Rv' 



and so using the no-shell crossing condition R' + Rv' . this 
is positive if and only if the numerator is positive. From 
(|45| . we have 



R't + Rtv' > R't 



Rt 



(recall that Rt < 0). The latter term is positive if and 
only if 



1 iR 

2 \~ 



-3/2 



1 



3"f 



1 



> 



F^ 
3F 



0(1). 



We note that 



1 



3"f 



111 
3/0 



0(r) > 0. 
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For t > tf^, we have 



-1 



Wc note that the coefficient of r here is positive. Hence 
if we can choose k so that 



tciO) /: 



6 /( 



i(K-A)-i = 



A 



3(A- k) 



>2, 



then there exists 5 > so that at > on k] as 
required. This choice entails k > |A, which can always 
be made. H. 

Proposition 6 Let k and 6 have the values required by 
Lemma 2. Then any future pointing polar null geodesic 
that enters the region 57 [(5, k] extends back to the central 
singularity in finite affine parameter time. 

Proof: It only remains to check the statement regarding 
finiteness of the value of the affine parameter s at which 
the geodesic meets the singularity. By Lemma 2 and 
(|4T|) . we see that i: < in the limit as the singularity 
is approached. If the limit t — > tc{0) is reached only as 
s — > — oo, then we would have 

lim i>0, 

s — > — oc 

contradicting the statement above. I 



VI. DISCUSSION 

We have revisited the issue of the visibility of the shell- 
focusing singularity in quasi-spherical dust collapse, mo- 
tivated by the observation that previous results have in- 
correctly assumed that there exist radial null geodesies 
in such space-times. As we have seen, this is not gen- 
erally the case: the existence of such geodesies implies 
an additional symmetry of the space-time. It is worth 
noting that our discussion has been restricted to the 4- 
dimensional case, whereas there have been several studies 
carried out in higher n + 2 = D > b dimensional Szek- 
eres space-times. We suspect that an analogous result 
applies: the existence of a radial geodesic in the space- 
time implies the existence of (n— 1) Killing fields, leaving 
just 2 + 1 non-ignorable coordinates. This conjecture is 



based on the structure of the derivates 



dv 
dxi ■■ 



1 < i < n in 



the higher dimensional case. However, we have not been 
able to determine the structure of the group of transfor- 
mations for the higher dimensional case that corresponds 
to the transformations ([8]) that play a crucial role in the 
proof of Proposition 1. 

It is perhaps worth pointing out that with the bene- 
ffi of hindsight, it is not surprising to see a connection 
between the existence of radial geodesies and symmetry. 
First, and on general grounds, it would be unusual to see 



a situation in which one had a conserved quantity (the 
values of the angles in this case) without the presence of 
some form of symmetry. Secondly and more specifically 
for this quasi-spherical situation, it is hard to envisage 
a geodesic emerging orthogonally from one 2-sphere and 
remaining orthogonal to other 2-spheres that it meets - 
unless the centres of those 2-spheres are aligned, with 
the alignment direction forming an axis of symmetry of 
the spacetime. This is exactly what we see happening in 
Proposition 1. 

In the axially symmetric case, there is one direction 
along which radial geodesies exist: the polar direction. 
We have looked briefly at the issue of the visibility of the 
singularity along this direction in the marginally bound 
case, and find no difference between the present case and 
the spherically symmetric case. However, we cannot gen- 
eralize our finding to either non-radial geodesies in the 
axially symmetric case, or to geodesies in the general 
case. The fundamental difficulty in doing so is that one 
cannot project the geodesic onto the r-t plane and retain 
all information required. The geodesic equations in the 
general case (no symmetry and hence no radial assump- 
tion allowed) form a second order nonlinear dynamical 
system with singular coefficients. The presence of the 
singular coefficients - which correspond to the points of 
space-time we are interested in analyzing - mean that 
standard methods of smooth dynamical systems do not 
offer means of approaching the problem. One possible ap- 
proach is to rescale the dynamical system by multiplying 
through by the most strongly vanishing denominator in 
the singular coefficients. One then absorbs this coefficient 
into a rescaled affine parameter, to obtain a smooth sys- 
tem. When this is done carefully, it is possible to convert 
the singular point of the original system to a stationary 
point of the rescaled system. However, this procedure 
typically yields non-hyperbolic equilibrium points and 
spurious equilibrium sets requiring the use of centre man- 
ifold analysis. Nonetheless, it has yielded useful results in 
a different context where similar problems (singular dy- 
namical systems) arise [l^ . It would be of interest to see 
if this approach could be used to study general geodesies 
in the non-axially symmetric Szekeres space-time, or in- 
deed the non-radial (non-polar) geodesies of the axially 
symmetric Szekeres space-time. The existence of addi- 
tional bounds on the metric functions (like (08])) would 
be of considerable use in this case. 
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